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statement  of  need 


Modern  rockets  operate  at  supercritical  pressures  with 
respect  to  the  propellants 

*  Need  to  understand  mixing  and  combustion  beyond 
liquid,  gas  states 

Shear  coaxial  injectors  are  a  common  choice  for  cryogenic 
liquid  rocket  engines 

Interactions  of  transverse  acoustics  with  injector’s  own 
modes  and  mixing  needs  to  be  understood  for  combustion 
instability 

*  Need  to  understand  differences  in  response  to  pressure 
and  velocity  nodes 

Understand  what  non-dimensional  numbers  capture  the 
mixing  of  typical  injectors 

Characterize  how  geometry  affects  mixing 
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Relevant  physics  of  shear  coaxial  jets 


1.  Transverse 
Acoustic  mode  from 
chamber/siren 

-  f=f(c,  geometry) 

2.  Acoustic  modes  for 
outer  and  inner  jets 

-  f~c/2L  -  2  speeds  of 
sound  up  and 
downstream 

3.  Wake  resuiting  from 
inner  post  thickness 

-  St=fw/U,, 

4.  Shear  iayer 
instabilities 

-  Ste=f0/U,, 

5.  Jet  preferred  modes 

-  St=fDij/Uij 


l_jp  Outer  Jet 
thickness 


Initial  Intermediate 

Merging  Zone  Merging  Zone 


Fuliy  Merge( 
Zone 


coaxiai  Jet 


wake 


shear  iayer 


Dahm  et  al,  JFM,  Vol  241, 1992 
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Highlights  from  literature  review  (1  of  2) 


•  Michaike,  1964 


—  Linear  stability  theory  for  inviscid  instabiiity  of  a 
hyperboiic  tangent  veiocity  profiie 

*  Crow  and  Champagne,  1971 

—  Singie  jet  preferred  mode,  St£|=fd/U~0.3 

*  Koetal,1 976-1 989 

—  Some  of  eariiest  detaiied  description  of  near  fieid  mixing 
for  coaxiai  Jets 

*  Boldman  et  ai,  1975 

—  Experimental  and  theoretical  analysis  for  mixing  of  two 
air  streams  with  different  velocities  -  points  out  different 
vortex  interactions,  St|~0.2  (Uave) 

*  Gutmark  and  Ho,  1983 

—  Collects  previous  results  on  Jet  preferred  mode,  St^  has  a 
range  from  ~0.24-0.64 


Kwan  and  Ko,  J.  Sound  and  Vibration,  48  (2),  1976 


Wicker  and  Eaton,  AiAA  J,  (32)  No.3,1994 


*  Wicker  and  Eaton,  1994 

—  Forces  air  inner  and  outer  jets  independently  -  observes 
vortex  growth 

*  Dahm  et  al,  1992 

—  Seminal  pictures  of  different  instabiiities  plus  effect  of 
absolute  velocity  and  R 
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Highlights  from  literature  review  (2  of  2) 


Balarac,  da  Silva,  Metals  et  al  (2003,  2007) 

—  DNS  analysis  of  coaxial  jets  -  same  density,  top-hat  profiles 

—  Consider  two  shear  layers,  study  effect  of  R 

—  Consider  axisymmetric  and  azimuthal  excitation 
Buresti,  Talamelli,  Petagna  (1994, 1998) 

—  Air  Jets,  same  density,  top-hat  profile,  Stj,o=fdo/Uoj~0.3  to  1, 
function  of  x/Dj 

Segalini,  Talamelli,  et  al  (2006,  2011) 

—  Air  Jets,  same  density,  top-hat  profile,  Stb(iip)=fb(lip)/U average 
Birbaud,  Ducruix,  Durox,  Candel  (2006-2007) 

—  Single  air  Jets,  low  Re,  laminar,  top-hat  profile,  subjected  to 
acoustic  modulation 

—  Systematic  study  of  effect  of  modulation  in  terms  of  Stj,,  Stg 
Tshohas,  Canino,  Heister  (2004,  2009) 

—  2D  unsteady  CFD  for  LOX/H2  elements  but  non-reacting 

—  Unforced  behavior,  found  St|ip=fd|ox/U|ox~0.1 0-0.25 
Richecoeur,  Scouflaire,  Ducruix,  Candel  (2006) 

—  Forced  transverse  acoustic  excitation  of  flames 


TjW 


L. 


Balarac  et  al,  Phys  of  Fluids  (19),  2007 


Segalini  et 
al,  Phys  of 
Fluids  (23), 
2011 


1 


Birbaud  et 
al,  Phys  of 
Fluids  (19), 
2007 


Richecoeur  et  al,  JPP  (22)  No  4,  2006 
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metric  Rationale 


LAR_ThinLip 


SAR_ThinLip 


u”  i'  .'i 


Two  recesses 


LAR_thickLip 

SAR_thinLip 

SAR_thickLip 

LARJhinLip 


D1 

(mm) 

D2 

(mm) 

D3 

(mm) 

D4 

(mm) 

t/D 

t={D2-D1)/2 

Ao/Ai 

0.51 

1.59 

2.42 

3.18 

1.05 

12.9 

1.40 

1.65 

2.44 

3.94 

0.09 

1.6 

1.47 

3.96 

4.70 

6.35 

0.84 

2.9 

0.70 

0.89 

2.44 

3.94 

0.13 

10.6 

SAR_ThickLip 


SAR,  LAR  ->  Small,  Large  Area  Ratio 
ThickLip,  ThinLip  ->  Post  lip  thickness 
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other  parameters  varied 


Momentum  flux  ratio 
Velocity  ratio 
Thermodynamic  state 
—  Subcritical  to  supercritical 


J=PoUo2/p|U|2 

R=Uo/U| 


Acoustics 


—  Acoustics  off 

—  Acoustics  on  (with  and  without  mode  matching) 

*  Phase  of  the  acoustics: 

—  0  degrees:  Pressure  Anti  Node  (PAN)  (velocity  node  (VN)) 

—  180  degrees:  Pressure  node  (PN)  (velocity  antinode  (VAN)) 


0 

PAN 

VN 


180 

PN 

VAN 
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Chronological  progression  (only  coaxial  results 
are  summarized  in  what  follows) 

^ - 

*  Single  jets,  no  coaxial  flow 

—  Davis  et.  al.  (Ph.D.  thesis)  -  singie  jets,  no  coaxiai  fiow 

*  Not  briefed  today 

*  Coaxial  Jets 

—  Leyva  et.ai.  -  LAR_thick 

—  Rodriguez  et.  ai.  (Ph.D.  thesis)  -  LAR_thick,  SAR_thin 

—  Graham  et.  al.  -  SAR_thin,  two  recesses 

—  Teshome  et.  ai.  (Ph.D.  thesis,  expected  Mar  2012)  - 
complete  all  four  geometries 

*  Aiso  complete  modal  analysis  of  eariier  geometries 

*  Future:  Combusting  coaxial  jets 

—  Wegener  et.  al.  (Ph.D.  thesis)  (in  process) 

—  Article  219  facility  funds 
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Image  interpretation  key 


Baseline: 

Acoustics 

off 


0“ 

Acoustics  on 
PAN,  VN 

“Pressure  coupled” 


180° 

PN,  VAN  - 
Largest  difference 

O 

expected  from  0 
“Velocity  coupled” 


pressure  =  fixed 
J  (poUo2/piUi2)  =  fixed 

PN  -  pressure  node  -  Min 
PAN  -  pressure  antinode  -  Max 
VN  -  velocity  node 
VAN  -  velocity  antinode 
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Baseline 


PAN 

p 

antinode 


Sub-Critical  Pressure:  Two  Geometries 

LAR_thickLip  LARJhinLip 


V- 


PN 

Pnode 


J=Po'Jo^/PiU|^ 

R=Uo/U| 


J  =0.2 


/?  =  2.0 
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Baseline 


PAN 


PN 


Near-Critical  Pressure:  Two  geometries 

LAR thickLip_ LARJhinLip 
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SAR_thinLip 


•  Nearcritical  conditions 

*  No  jet  bending 
observed 
•  large  vortical 
structures 
generated  when 
coupled  to  injector 
mode 

*  Reduction  of  dark 
core  can  be  as  iarge 
as  90% 

*  More  ciear  response 
of  Jet  to  pressure 
antinode 

*  Subcriticai  conditions 

□  •  Same  mode  - 

vorticai  structures 
*  Not  as  dramatic 
reduction  as  with 
nearcriticai  cases 


(~)  Pr=1.05,  J=0.5,  R=2 
^^Pr=1.05,  J=0.5,  R=1.3 
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Pr=1.05,  J=1.7,  p7p=0.32% 
Coupled  to  Outer  Jet  Acoustic  Mode 


Pr=1.05,  J=2.2, 
p7p=0.25% 

Non  Coupled 


Effect  of  Recess:  SAR_ThinLip;  Pr=0.45,  J~0.09 


:at 


Flush  inner  post,  J=0.09,  p7p=0.45% 


Recessed  inner  post,  J=0.089,  p7p=0.60% 


Qualitatively  similar  at  very  low  J  values 
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^  Acoustic  Analysis  for  Injectors 


In  collaboration 
with  Jeff  Muss  and 
Rory  Davis,  Sierra 
Engineering 


Have  more  accurately 
computed  the  acoustic 
modes  for  the  inner  and 
outer  jets  for  constant  and 
linearly  varying 
temperatures  for 
subcritical  and 
supercritical  pressures 
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Amplitude  Amplitude  Amplitude 


New  analysis:  synchronized  p’  and 
images  taken  with  microscopic  iens 


1  p 

0.9  - 
0.8  - 
0.7  - 
0.6  - 
0.5  - 
0.4  - 
0.3  - 
0.2  - 
0.1  - 
0  — 


Chamber 
Center 


3000  6000  9000  12000 

Frequency  (Hz) 


Inner  Jet  Pressure 


Chamber  Wall  Center  Pressure 


n. 

"S. 


Outer  Jet  Pressure 


Pr  =  1.03,  J  =  5,1,VR  =  5.9 
Forcing  Condition:  3000  Kz,  Maximum  Ap 


3000  6000  9000  12000 

Frequency  {Hz} 


Uo=6.  25  m/s,  Ui=1.06  m/s,  To=165K  ,  Ti=118  K 

Distribution  A:  Approved  for  public  release;  distribution  unlimited 


Case  I:  Baseline  Flow 


Snapshot  Average 


Mode  1  Mode  2  Mode  3  Mode  4 


Mode  5 


Case 

Pr 

^R,OJ 

^R,IJ 

R 

J 

I 

0.44 

1.19 

0.84 

5.8 

2.0 

Mode  6 


i 
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Case  I:  Acoustically  Forced  Flow 


Mode  1  Mode  2  Mode  3 


Mode  4 


Mode  5  Mode  6 


Case 

Pr 

^R,OJ 

^R,IJ 

R 

J 

/(kHz) 

I 

0.44 

1.19 

0.84 

5.8 

2.0 

2.96 
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Collaboration  with  ECP:  Thomas  Schmitt,  Juan  Rodriguez 
(AFRL  post-doc),  Sebastien  Candel,  ivett  Leyva  -  enirau 

19 

'  Velocity  Fields  for  Different  J  values  ” 


600  = 


500  - 


.^300-  t 


iOO-  ’ _ _ 


Longitudinal  slice  of  axial  velocity.  White:  minimum;  black:  maximum. 
Dark  line  indicates  iso-contour  of  zero  axial  velocity.  Top:  J  =  1.1;  middle:  J 

=  3.0;  bottom:  J  =  9.3. 


Top:  centerline  profile  of  density;  bottom: 
centerline  profile  of  velocity.  Dark  line:  J  = 
1.1;  dash  point  line:  J  =  3.0;  dashed  line:  J  = 
9.3. 
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^  Fundamental  frequencies  for 
baseline  conditions 


Case  “N6” 
(J=3.05) 
Density 
distribution 
(white:  60 
kg/m^; 
biack:  410 
kg  /m^; 
iogarithmic 
scaie). 


OUTER 

JET 

J=1.05 

R=2.48 


+ - * 


S  1  t  £  S  E  E  \ 

xjfl  M 

^ai: 


Probe  A 


Fixqunin-  (Hi) 


MB?  r:iMO 

Frvqttnr;;'  Ufa) 


Case 

St«=hef/U,j 
(probe  A) 

St*  =dif/Uij 
(probe  B) 

St'=l|f/U.,j 
(probe  C) 

N2 

0.25 

0.34 

0.15 

N6 

0.25 

0.26 

0.14 

Case  “N8” 

(J  =  9.3) 
Density 
distribution 
(white:  60  kg/ 
m^;  biack: 
maximum ; 
iogarithmic 
scaie). 


OUTER 

JET 

J=3.05 

R=4.18 


Found  relevant  St  numbers  for  our  configurations 
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Conclusions 


For  LARthickLip 


—  Found  bending  mode  with  iargest  effect  at  velocity  antinodes 

—  The  reduction  on  the  dark  core  length  was  greatest  for  a  medium  J  range 

—  For  near  critical  pressures,  the  collaboration  with  ECP  determined  relevant  St  for  our 
injector  configuration  and  was  able  to  capture  qualitative  behavior  of  natural  and  excited 
jets 

For  SARthinLip 

—  Did  not  see  bending  mode  for  conditions  studied 

—  Saw  vortex  roll-up  and  puffing  occurring  over  entire  J  range  (0.09-21)  tested  -  most 
predominantly  at  pressure  nodes  and  most  likely  associated  with  acoustic  mode  of  outer 
jet 

—  Vortex  roll  up  and  puffing  is  independent  of  injector  recess  configuration 

For  LARthinLip 

—  Sees  both  bending  and  vortex  roll-up  modes  depending  on  the  acoustic  frequency 


Synchronized  measurements  for  p’  and  jet  images  -  correlate  structures 
observed  to  p’  amplitudes 


Initial  POD  analysis  recovers  features  seen  in  movies  -  next  wiii  be  to 
compute  convective  veiocities  of  structures  observed 
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AFOSR/NASA  Combustion 
Workshop 


•  The  single  largest  unknown  in  combustion  instabiiity  is  the 

“combustion  response”  (how  combustion  responds  to  acoustic  waves) 


•  Within  the  combustion  response,  flame  holding  in  the  near  injector  fieid 
is  a  key  mechanism 


•  AFRL  -  Edwards  was  tasked  to  deveiop  an  experiment  and  to  lead  a 
joint  experimental  /  modeling  team  to  study  coupled  flame  holding 
mechanisms. 


-  LEVERAGE  EXPERIENCE  GAINED  WITH  COLD  FLOW  EXPERIMENTS 
•  GA  Tech  was  seiected  to  iead  an  effort  on  a  ciosed-loop  study. 


Stability 
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EC-4H  -  Combustion  Instability  La 


Concentrate  on  near  injector  field 


Measure  p’  and  q’  simultaneously  to 
evaluate  the  Rayleigh  criteria  for 
combustion  instability 


Start  with  shear  coaxial  jets  -  cold  flow 
heritage 


Lab  designed  for  2000  psi  -  about 
double  the  pressure  from  other  labs  in 
the  world 


Start  with  current  design  for  acoustic 
drivers  pressure  nodes  and  antinodes 
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BACKUP 
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Image  interpretation  key 


pressure  =  fixed 


J 


(PoUo2/p.Ui2)  =  fixed 


Baseline: 

Acoustics 

off 


180° 

PN,  VAN  - 
Largest  difference 

O 

expected  from  0 
“Velocity  coupled” 


Acoustics  on 


0° 

PAN,  VN 

“Pressure  coupled” 


Subsequent  images: 
acoustics  on; 

O 

Phase  increases  by  45 


PN  -  pressure  node  -  Min 
PAN  -  pressure  antinode  -  Max 
VN  -  velocity  node 
VAN  -  velocity  antinode 


360° 

should  look 
similar  to  0° 
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newl 


new3 


^  /  RnewJ 


4  A 

Geometries  for  present  study 


SAR_thinLip 


=  0,70  mm 
=  0,83  mm 
=  1,22  mm 
=  1,97  mm  ,■ 


R21 

R22 

R23 

R24 


A 

m  • 

9- 

^  *  ’ 

j 

V 

^  ' 

-ip 

I 


LAR=LargeAreaRatio 

SAR=SmallAreaRatio 


D1  D2  D3  D4  t/D  Ao/Ai 

(mm)  (mm)  (mm)  (mm)  t=(D2-Di)/2 


LAR_thickLip 

0,51 

1.59 

2.42 

3.18 

1.05 

12.9 

SAR_thinLip 

1,40 

1.65 

2.44 

3.94 

0.09 

1.6 

SAR_thickLip 

1,47 

3.96 

4.70 

6.35 

0.84 

2.9 

LAR_thinLip 

0,69 

0.89 

2.44 

3.94 

0.15 

ll.Oi 
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4  t 

Relevant  variables  for  cold-flow  studies 

^ 


Geometry 

Acoustics 

Recess 

Phase 

mm 

Coupling 

Single  Jet 

On/Off 

N/A 

v’  max 

N/A 

No 

LAR_thickLip 
(injector  1) 

1/2D1 

No 

SAR_thinLip 
(injector  2) 

Off 

P’  max 

U’max 

1/2D1 

2-phase  P<Pc 

^  0.1-20 

SAR_thickLip 
(injector  3) 

P>Pc  0.1-20 

^  T>Tc 

°  T<Tc 

Yes&  No 

LAR_thinLip 
(injector  4) 

0 
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